statics and dynamics of elastic manifolds in media with long-range correlated disorder 
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We study the statics and dynamics of an elastic manifold in a disordered medium with quenched defects 
correlated as ~ for large separation r. We derive the functional renormalization-group equations to one- 
loop order, which allow us to describe the universal properties of the system in equilibrium and at the depinning 
transition. Using a double e — 4 — d and 5 — 4 — a expansion, we compute the fixed points characterizing 
different universality classes and analyze their regions of stability. The long-range disorder-correlator remains 
analytic but generates short-range disorder whose correlator exhibits the usual cusp. The critical exponents and 
universal amplitudes are computed to first order in e and S at the fixed points. At depinning, a velocity-versus- 
force exponent /3 larger than unity can occur We discuss possible realizations using extended defects. 



I. INTRODUCTION 

Elastic objects in random media are the simplest example 
of a disordered system exhibiting metastability, glassy behav- 
ior, and dimensional reduction, which are difficulties present 
in a broader class of disordered systems d |^ [H. They can 
be used to model a remarkable set of experimental systems. 
Domain walls in magnets behave as elastic interfaces and can 
experience either random bond disorder (RB) as in ferromag- 
nets with nonmagnetic impurities, or random field disorder 
(RF) as in disordered antiferromagnets in an external mag- 
netic field |4]. The interface between two immiscible liquids 
in a porous medium exhibits the same behavior and under- 
goes a depinning transition as the pressure difference is in- 
creased liSj]. Charge-density waves (CDW) in solids show a 
similar conduction threshold J^. Another example of peri- 
odic systems are vortex lines in superconductors which can 
form different glass phases in the presence of weak disorder 
IItI HI nil- In all these systems, the interplay between elastic 
forces that tend to keep the system ordered, i.e., flat or pe- 
riodic, and quenched disorder, which promotes deformations 
of the local structure, forms a complicated energy landscape 
with numerous metastable states. This results in glassy prop- 
erties and a nontrivial response of the system to external per- 
turbations. In particular, the interface becomes rough with 
displacements growing with the distance x as 



(1) 



where ( is the roughness exponent. Elastic periodic structures 
in the presence of disorder lose their strict translational or- 
der and form quasi-long-range order characterized by a slow 
growth of displacements, 



C{x) = Adl-o-lxl 



(2) 



where the amplitude Ad is universal in the simplest case. At 
zero temperature, a driving force / exceeding the threshold 
value fc is required to set the elastic manifold into steady mo- 
tion with a velocity v that vanishes as v ^ (/ ^ fc)^ at the 
transition point. The correlation length diverges close to the 
transition f — fc £, ^ {f — fc)^" and the characteris- 
tic time as T ^ where z is the dynamic critical exponent. 
Note that the roughness exponent and the universal amplitudes 



determined at the depinning transition are in general different 
from the exponent and amplitudes measured in equilibrium. 

Two methods were developed to study the statics of an elas- 
tic manifold in a disordered medium. One of them is the Gaus- 
sian variational approximation (GVA) performed in replica 
space, which can be applied to both classes of elastic man- 
ifolds, i.e., to interfaces [10] without overhangs and to peri- 
odic systems (SdH]. Within this approach, which is believed 
to be exact in the mean-field limit, i.e., when the manifold 
lives in a space of infinite dimensions, metastability is de- 
scribed by breaking of replica symmetry, which allows one 
to compute the static correlation functions and to obtain dif- 
ferent thermodynamic properties. Another method that can 
be applied to dynamics as well as to statics is the functional 
renormalization group (ERG) |fl2ll . Simple scaling arguments 
show that large-scale properties of elastic systems are gov- 
erned by disorder for d < d^c — 4 and that perturbation 
theory in the disorder breaks down on scales larger than the 
so-called Larkin scale [13]. To overcome this difficulty, one 
performs a renormalization-group analysis. It was shown that 
in this case one has to renormalize the whole disorder corre- 
lator that becomes a nonanalytic function beyond the Larkin 
scale fl^ [13, [m [l^. The appearance of a non-analyticity 
in the form of a cusp at the origin is related to metastability, 
and nicely accounts for the generation of a threshold force at 
the depinning transition. It was recently shown that the ERG 
can unambiguously be extended to higher loop order so that 
the underlining nonanalytic field theory is probably renormal- 
izable to all orders [17. .18,, .19.]. Although the two methods, 
GVA and ERG, are very different, they provide a fairly con- 
sistent picture of the statics, and recently a relation between 
them was found [20]. There is also good agreement with 
results of numerical simulations, not only for critical expo- 
nents [ 2TII22II23I1 but also for the whole renormalized disorder 
correlator ||24|| . However, many questions remain open. Al- 
though the dynamics in the vicinity of the depinning transition 
and at zero temperature is well understood, there is no satis- 
factory theory for finite temperature, and in particular for the 
thermal rounding of the depinning transition [25]. It is also re- 
markable that the exponent f3 in experiments on depinning is 
usually larger than 1, while ERG and numerical simulations of 
elastic systems with weak disorder give values smaller than 1. 

Most studies of elastic manifolds in a disordered medium 
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treat uncorrelated pointlike disorder Real systems, however, 
often contain extended defects in the form of linear disloca- 
tions, planar grain boundaries, three-dimensional cavities, etc. 
It is known that such extended defects, or pointlike defects 
with sufficiently long-range correlations, can change the bulk 
critical behavior 0261 l27L Flux lines in su- 

perconductors are the most prominent example. The pinning 
of the flux lines by disorder prevents the dissipation of energy 
and determines the critical current Jc, which is of great impor- 
tance for applications. It was found that extended defects pro- 
duced, for instance, by heavy-ion irradiation, can increase 
by several orders of magnitude |33]. Systems with anisotropic 
orientation of extended defects can be described by a model 
in which all defects are strongly correlated in ed dimensions 
and randomly distributed over the remaining d ~ £d dimen- 
sions. The case = is associated with uncorrelated point- 
like defects, while extended columnar or planar defects are 
related to the cases £d = \ and 2, respectively. The bulk- 
critical behavior in the presence of this type of disorder was 
studied in Refs. ll26ll27il28i l29ll using a perturbative RG anal- 
ysis in conjunction with a double expansion in e = 4 — d and 
Ed- The pinning of flux lines by columnar disorder was stud- 
ied in Ref. i34il . where it was shown that the system forms 
a Bose glass phase with flux lines strongly localized on the 
columnar defects, resulting in a zero dc linear resistivity. It 
was argued recently that the topologically ordered glass phase 
(Bragg glass) formed by flux lines can be destroyed in the 
vicinity of a single planar defect ll35ll . It has been shown that 
the small dispersion in orientation of columnar defects forms a 
new thermodynamic phase called "splayed glass" 1 36] . In this 
phase, the entanglement of flux fines enhances significantly 
the transport of superconductors fyf\. Competition between 
various types of disorder, point and columnar, has also been 
studied, at equilibrium i38[[39tl and in the moving phases ll43l . 

In the case of an isotropic distribution of disorder, power- 
law correlations are the simplest example with the possibility 
for a scaling behavior with new fixed points (FPs) and new 
critical exponents. The bulk-critical behavior of systems in 
which defects are correlated according to a power law 
for large separation r was studied in Refs. ifsol [3TI l32ll . The 
power-law correlation of defects in d-dimensional space with 
exponent a — d — Ed can be ascribed to randomly distributed 
extended defects of internal dimension Ed with random orien- 
tation. For example, a ~ d corresponds to uncorrelated point- 
like defects, a = d— l(a = d— 2) describes infinite lines 
(planes) of defects with random orientation. In general, one 
would probably not expect a pure power-law decay of correla- 
tions. However, if the correlations of defects arise from differ- 
ent sources with a broad distribution of characteristic length 
scales, one can expect that the resulting correlations will 
over several decades be approximated by an effective power 
law fid]. If the correlation function of disorder can be ex- 
pressed as a finite sum of power-law contributions Qr""', 
one can expect that the scaling behavior is dominated by the 
term with the smallest [30]. Power-law correlations with a 
noninteger value a = d~ df can be found in systems contain- 
ing defects with fractal dimension df 1141(1 . For example, the 
behavior of ''He in aerogels is argued to be described by an XY 



model with LR correlated defects This is closely related 
to the behavior of nematic liquid crystals enclosed in a single 
pore of aerosil gel, which was recently studied in Ref. ll43ll . 
using the approximation in which the pore hull is considered 
a disconnected fractal. Finally, studies of the Kardar-Parisi- 
Zhang (KPZ) equation with power-law correlations in time 
bear connections to the case d ~ 1 considered here. 
However, the perturbative method used there cannot address 
directly the zero-temperature (strong KPZ coupling) phase, 
contrary to our present study. 

In the present paper, we study the statics and dynamics of 
elastic manifolds in the presence of (power-law) LR correlated 
disorder using the FRG approach to one-loop order. The paper 
is organized as follows. Sectionllll introduces the model. Pos- 
sible physical realizations are considered in Sec. |III] Section 
HVl describes the dynamical formalism and perturbation the- 
ory. In Sec. [V] we renormalize the theory and derive the FRG 
equations to one-loop order. In Sec. [VT] we study random 
bond, in Sec. IVIH random field, and in Sec. IVIIII periodic 
disorder In Sec. |IX]we discuss fully isotropic extended de- 
fects. In the final section, we summarize the obtained results 
and our conclusions. 



II. THE MODEL 

We consider a d-dimensional elastic manifold embedded 
in a D-dimensional space with quenched disorder. The con- 
figuration of the manifold is described by an A^-component 
displacement field denoted below u{x), or equivalently Ux, 
where x denotes the d-dimensional internal coordinate of 
the manifold. For example, a domain wall corresponds to 
d = D - 1 and iV = 1, a CDW to d = £> and = 1, 
and a flux lattice to d = D and = 2. In what follows, we 
focus for simplicity on the case = 1 and elastic objects with 
short-range elasticity. Extensions to A^ > 1 and LR elasticity 
are straightforward for the statics. The energy of the manifold 
in the presence of disorder is defined by the Hamiltonian 



d"x 



-[Wu{x)]^ + V{x,u{x)) 



(3) 



where c is the elasticity and y is a random potential. In this 
paper, we study the model where the second cumulant of the 
random potential has the form 



V{x,u)V{x',u') 



Ri{u-u')5'^{x-x') 
+R2{u-u')g{x-x'). (4) 



The first part corresponds to pointlike disorder with short- 
range (SR) correlations in internal space. The second part cor- 
responds to long-range (LR) disorder in internal space and the 
function g{x) ^ x^^ at large x with a > 0. For convenience, 
we normalize it so that its Fourier transform is g{q) = \q\°'~'^ 
at small q with unit amplitude. A priori we are interested in the 
case a < d, where the correlations decay sufficiently slowly 

in internal space. We denote everywhere below = J J27t)^ 
and — J d'^x. The short-scale uv cutoff is implied at q ^ A 
and the size of the system is L. 
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One could start with model (|4|l, setting i?i = 0; however, 
as we show below, a nonzero i?i is generated under coarse 
graining. Note that the functions Ri{u) can themselves a pri- 
ori be SR, LR, or periodic in the direction of the displacement 
field u. For SR disorder in internal space only, i.e., R2 — 0, 
these cases are usually referred to as random bond (RB), ran- 
dom field ^ \u\ at large m] (RF), and random periodic 
(RP) universality classes. Below we discuss how these classes 
extend to the case of LR internal disorder {R2 nonzero). 

The model (|3]l and (HJl could easily be studied using 
presently available numerical algorithms for directed mani- 
folds, in its statics (e.g., exact ground-state determinations) 
and its dynamics (e.g., critical configuration at depinning), by 
directly implementing a random potential correlated as de- 
scribed by Eq. (HJi. It is also interesting to examine which 
type of correlations in a random medium can naturally lead to 
Eq. dUi and how such disorder could be realized from, e.g., 
distributions of extended defects, since some of them may be 
experimentally feasible. 

III. REALIZATIONS AND UNIVERSALITY CLASSES 




FIG. 1: (Color online) Linear defects randomly and isotropically dis- 
tributed on parallel planes with random distances between them. This 
geometry mimics distribution (|6). 

Although we mostly discuss extended defects, other 
sources of long-range correlations are possible, such as de- 
fects where each single one creates a long-ranged disorder 
potential, or a substrate matrix itself quenched at a critical 
point. 



A. Defect potential 

Let us first recall how long-range correlations can arise in 
the potential created by defects. To this purpose, call v{r) the 
defect potential, in the simplest case taken to be proportional 
to defect density. Consider for simplicity a large number of 
weak defect lines with a uniform and isotropic distribution in 
a space of dimension D. These create an almost Gaussian 
random potential v{r) with 



v{r)v{r') - r — for r -> 00 (5) 

and a = D — \. To derive this, consider defects of finite 
radius ad- The probability that point r' is contained in the 
defect going through r is ~ (orf/k ^ '"'D^^^. i-e-^ inversely 
proportional to the sphere of radius \r — r'\. This is easily 
generalized to isotropic distributions of extended defects of 
internal dimension e^;, with a — D~ed- Note that by extended 
defects we mean defects that are perfectly correlated along 
their internal dimension. Generalizations where defects are 
themselves (anisotropic) fractals can also be considered. 

An important case is a uniform distribution of extended de- 
fects in IJ-dimensional space, but isotropic only within a lin- 
ear subspace of dimension D' . For instance, one can irradiate 
a material in the bulk while simultaneously rotating it along an 
axis. This produces a distribution of linear defects {Sd — 1), 
isotropic within the plane (£)' = 2), and normal to the axis 
(see Fig. [Til. More generally, this yields a defect potential with 
second cumulant, 

w(r, z)v(r' , z) = g{r — r')f{z — z'\ (6) 
g(r) ~ r-'^ , 

while j{z) is short-ranged (here r e , z €E R^'^ , a = 
D' - Ed). 



B. Coupling to the manifold 

We now examine how the long-range correlated defect po- 
tential couples to the elastic manifold and what type of LR 
model results. A general formulation of this coupling (see, 
e-g-, yj) has the form 

v-fe (7) 

where the defect potential lives in the iJ-dimensional space 
parametrized by (x, u) and x E R'^ is the internal coordinate 
of the manifold. p{x, z, u) is the manifold density. Each type 
of coupling to the disorder corresponds to a different function 
p{x, z, u), and we now indicate the main cases. 

1. Elastic interfaces in random bond disorder 

Let US first discuss elastic interfaces in the so-called random 
bond (RB) case, where the coupling between disorder and in- 
terface occurs only in the vicinity of the interface as, e.g., for 
domain walls in magnets with random bond disorder This 
corresponds to the choice 

p{x, z,u) ^ 5{z — u), (8) 

hence the additional variable z introduced in Eq. ^ is identi- 
cal to u, the displacement field (with in general D — d = N ). 
In that case, 

Vrb{x,u) v{x,u). (9) 

Consider now a uniform distribution of defects in the D- 
dimensional plane but isotropically distributed within the ( av- 
eraged direction) of the internal space of the manifold. This 
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is given by Eq. ^ above with D' = d, 



Vrb{x,u)Vrb{0,0) ^ g{x)R2{u), (10) 

which is model ^ with a SR function i?2 (u) and, in full gen- 
erality, a = d — £d- The physical realization in terms of ex- 
tended defects is thus an interface (c? = 2) in D = 3 with line 
defects all orthogonal to the u directions, isotropically dis- 
tributed within the (average) plane of the interface, and a—1. 
This is illustrated in Fig. [T] 

Another physical realization consists of extended defects 
with finite random lengths such that the distribution of lengths 
has a power-law tail for large lengths. For instance, needles 
of variable lengths aligned along one direction could act on a 
directed polymer d = 1 as power-law correlated disorder in 
internal space. 

An interesting, though qualitatively different, case occurs 
when the extended defects are distributed isotropically in the 
whole (x, u) space. This will be discussed in Sec. |lXl Finally, 
note that we consider weak Gaussian disorder It is possible 
that at strong disorder another phase exists where the line or 
manifold gets localized along the strongest extended defect. 



2. Elastic interfaces in random field disorder 

Random field (RF) disorder is described by the function 

p{x, z,u) ^ Q{u — z), (11) 

where Q{z) is the Heaviside step function. This means that 
the change in energy when the interface moves between two 
configurations is proportional to the sum of all defect poten- 
tials in the volume (in R^) spanned by this change. The dis- 
cussion of the geometry of defects needed to produce LR dis- 
order in internal space is identical to the last section. Substitu- 
tion of Eq. ( fTTT ) into Eq. (|7]) yields the RF disorder correlator, 
which can be approximated by Eq. (|4|i with Ri{u) ^ — u for 
large u. 



3. Periodic systems 

As an example of periodic systems, we consider incommen- 
surate single-2 CDWs. In that case D ^ d, hence the func- 
tion p{x, z, u) = p{x, u) in Eq. (|7|. The electron density of 
CDWs neglecting effects caused by an applied strain has the 
form [3, 6] 

p{x, (/)) ~ po + Pi cos{2A;F[a::j_ - "(a;)]}, (12) 

where the displacement u{x) of the maximum of the density 
is related to the standard phase field via <j){x) — —2kpu{x), 
where fcp is the Fermi wave vector. The d-dimensional space 
is split into x — , x±), with x± denoting the modulation 
direction of the CDW and fcp the Fermi wave vector. 

We again consider the situation of extended defects all 
aligned with the direction x\^ and isotropically distributed in 



that subspace. The random potential experienced by the CDW 
is given by 

V{x,(j)) = hi{x) cos (j){x) + h2{x) sm(t){x), (13) 

with Gaussian distributed hi{x) — v{x) cos{2kpx±) and 
h2{x) — v{x) sm{2kpx±). On large scales kFX± ^ 1, and 
their cumulant can be approximated by [from Eq. (|6])] 

'RxWm - 14rS^,S\x) + ^^S.,S{x^), (14) 

where we have omitted all rapidly fluctuating contributions. 
Equations (fTST l and ( fT4] i give the potential correlator in a form 
that can be generalized to 

V{x,u)V{x',u') = Ri{u- u')S'^{x - x') 

+R2iu ~ u')g{xi\ -x\)S'^^{xj_-x'^), (15) 

with d± ^ 1 and bare functions Ri{(j>) = ^vf cos(0), u = (p. 
Thus periodic systems are described by periodic functions 
Ri{u). Here is the dimension of the transverse subspace. 
Note that the Hamiltonian Hxy = J d'^x[^{'V(f))'^ + V{x, 0)] 
with V{x, (f) given by Eq. ( fT3] l and a Gaussian distribution of 
fields hi{x)hj{x') ~ g{x — x') describes the XY model with 
long-range correlated random fields. Therefore, the latter can 
be mapped onto periodic manifolds with correlator (flSl l and 
d± = 0, i.e., to model (|4|i with periodic functions Ri{u). In 
the next section, we will show how the ERG picture of model 
( fTSb can be obtained from the ERG results for model (|4|i. It 
is worthwhile to note that in the case of periodic systems, the 
integration in Fourier space is supposed to be over the first 
Brillouin zone. Note also that we have neglected the cou- 
pling of disorder to the long wavelength part of the density 
— Po / d'''xv{x)Wu{x) as it is usually irrelevant near the upper 
critical dimension. Indeed, in the replicated Hamiltonian (see 
below), this coupling generates additionally to the SR term 
-1/T J d'^xai\'uaix)\'ui,{x) the LR term 

J d'^xd''x'a2gixii-x[)S''^{x±-x'^)'^Ua{x)\'ub{x'). 

For small disorder in the vicinity of the upper critical dimen- 
sion, both of them renormalize to zero according to 

deai = (2-rf-2C)ai + ..., (16) 
de(J2 = (2-a-d_L-2C)o-2 + ... (17) 

IV. DYNAMICAL FORMALISM 

The overdamped dynamics of the elastic manifold in a dis- 
ordered medium can be described by the equation of motion 

rjdtUxt = cV^Uxt + F{x, Uxt) + fxt, (18) 

where -q is the friction coefficient. In the presence of 
an applied force /, the center-of-mass velocity is w = 

L^"^ J^dtu^f The pinning force reads F = —duV{x,u), 



and thus, for correlator (|4|i, the second cumulant of the force 
is given by 



Ai{u-u')5''{x-x') 
+A2iu-u')gix~x'), (19) 



F{x,u)F{x',u') 



with Ai = —R'/{u) in the bare model. In the following, we 
will always use g{q) = and g{x) = e'"^^g{q). 

The most important quantity of interest is the roughness ex- 
ponent ( measured in equilibrium or at the depinning transi- 
tion f = fc defined by 



C{x~x') = \i 



r'|2C 



(20) 



t. — i 

b 
c 




FIG. 2: Diagrammatic rules: a, propogator; b, SR disorder vertex; 
c, LR disorder vertex; d and e, one-loop diagrams generating the 
critical force at the depinning and giving correction to the mobility 
and elasticity. 



The velocity vanishes at the depinning transition as v ^ \ f — 
fc\^, while the correlation length diverges at the transition as 
C ~ 1/ ~ fc\~'^- One can also introduce the dynamic critical 
exponent z, which relates spatial and temporal correlations via 

t - a;^ 

Let us briefly sketch how one can construct the perturba- 
tion theory in disorder We adopt the dynamic formalism. It 
also allows us to obtain the statics equations (to one loop and 

= 1 these can easily be deduced, as can be checked using 
replica). Instead of a direct solution of the equation of motion 
( fTSl l with consequent averaging over different initial condi- 
tions and disorder configurations, we employ the formalism 
of generating functional. Introducing the response field iixt 
we derive the effective action, which reads 



S = 



iuxtivdt - cV + m )uxt - / iuxtfxt 

xt J xt 

i / iuxtiuxt' Ai{uxt - Uxt') 

^ JxtV 

^ iuxtiux't'g{x - x')A2{uxt- Ux't')A'2-l) 



where we have added a small mass m, which plays the role 
of an IR cutoff. To study the critical domain, one has to take 
the limit m ^ 0. The average of the observable over 
dynamic trajectories with different initial conditions and over 
different disorder configurations can be written as follows: 



{A[uxt]) = / V[u]V[u]A[uxt]e- 



-S[uM] 



(22) 



Furthermore, the response to the external perturbation fxt, 
which is local in time and in space, can be computed using 
{A[uxt]iuxt) — jj-^ i^b^xt]) ■ Note that causality is fulfilled, 
and here we adopt the Ito convention, which results in getting 
rid of all closed loops composed of response functions. 

In the absence of LR correlated disorder action, Eq. (ISTT i 
exhibits the so-called statistical tilt symmetry (STS), i.e., the 
invariance of the disorder terms under the tilt Uxt Uxt + hx 
with an arbitrary function hx- The STS gives the exact identity 
J^Tlqt = ^/cq^ for the response function 7?.qt — {uqt iu-qo), 
which implies that the elasticity is uncorrected by disorder to 
all orders. LR correlated disorder destroys the STS of action 
tni . and thus, in principle, allows for a renormalization of the 



elasticity. The quadratic part of the action ( |2TI ) yields the free 
response function 



[Uqt m-go)o = Rqt = 



7] 



{cq^+rn'^)t/ri 



(23) 



which can be used to generate the perturbation theory in dis- 
order. The theory has two disorder interaction vertices Ai{u) 
and A2{u). At each vertex Ai{u) there is one conservation 
rule for momentum and two for frequency while each vertex 
A2 (m) carries additional momentum dependence. In what fol- 
lows, we generalize the splitted diagrammatic method devel- 
oped in Ref. fis'], shown in Fig. |2l As is the case for the 
model with SR disorder, our model exhibits the so-called di- 
mension reduction, both in the statics and in the dynamics. 
The naive perturbation theory obtained taking the functions 
Ai{u) analytic at u = leads to the same result as that com- 
puted from the Gaussian theory setting Ai{u) — Ai{Q). In the 
limit m ^ 0, the two-point function then reads to aU orders: 



UgtU-qt' 



Ai(0) , A2(0) 



(24) 



The first term in Eq. (l24l i dominates in the limit g — > for 
a > d, and LR disorder is irrelevant in this case, while the 
last term dominates for a < d. Equation (l24l l results in 
^ = (4 - d)/2 for a > d and C = (4 - a)/2 for a < d, which 
are known to be incorrect. The physical reason for this is the 
existence of a large number of metastable states. The rough- 
ness exponent can be estimated using Flory arguments setting 



Then the gradient term scales as V 



The 



pinning force for SR disorder scales as F{x, Ux) ^ 2:"^'*+^^/^ 
and for LR disorder as F{x, Ux) ^ Therefore, in 

the regime where the behavior is governed by SR disorder, 
the Flory estimate gives for RF disorder the Imry-Ma value 
CsR ~ (4 ^ c?) /3 while for LR RF disorder we get = 
(4 — a)/3. A similar argument constructed from the potential 
correlators Ri{u) yields the Flory estimates ^gj^ = (4 — d) /5 
and Clr — (4 ~ a)/5, respectively, for the case of random 
bond disorder To obtain corrections to the Flory values, the 
ERG developed in Refs. d H [3 [M HI, [3 

will be em- 
ployed. The solution is nontrivial because the renormalized 
disorder becomes nonanalytic above the Larkin scale, and one 
has to deal with a nonanalytic field theory. Here we generalize 
this approach to the case of LR correlated disorder. 



V. FUNCTIONAL RENORMALIZATION 

We now consider the renormalization of model (1211 1. The 
subtleties arising for the correlator ( fTsT i will be discussed 
briefly at the end. We carry out perturbation theory in the 
bare disorder correlators A^o (u) and then introduce the renor- 
malized correlators Ai{u). We will suppress the subscript "0" 
to avoid an overly complicated notation. According to the 
standard renormalization program, we compute the effective 
action to one-loop order. Here we adopt the dimensional reg- 
ularization of integrals and employ the minimal subtraction 
scheme to compute the renormalized quantities and absorb the 
poles in £ = 4 — d and 6 = A — a into multiplicative Z factors. 
When derivatives of the at m = occur, in the dynamics 
(i.e., at the depinning transition for dynamical quantities) they 
are taken at u = 0+ as can be justified exactly for N — 1. In 
the statics, the treatment is more subtle (as discussed in two- 
loop studies lO) but is not needed in the present one-loop 
study. 

Let us firstly consider the first-order terms generated by ex- 
pansion of in disorder. These terms are given by diagrams 
d and e shown in Fig.|2] We start from 



t>t',x 



iUxtA2{uxt - Ux't')g{x - x')iUx't'-(25) 



t>t',x,x' 



Expanding A^ (u) in a Taylor series and contracting one m, we 
obtain the leading corrections to the threshold force, friction 
and elasticity. The terms giving the threshold force to leading 
order are 



iUxt^iiO )Rx=o,t-t' 



t>t',x 



iu^tA2(0+)5(a; - x')Rx~x' ,t~t' ■ (26) 

t>t' ,x,x' 

They are strongly uv diverging (~ A''"^ + A°~^), and thus 
are nonuniversal. The terms proportional to A"(0+) can be 
rewritten as corrections to friction and elasticity using the ex- 
pansion 

Uxt - Ux't' = {t- t')dtUxt + {x- x')i——Uxt 

OXi 

+ (x - x'Ux - x'),\p^ + 0{M\Ax^). (27) 
The first term in Eq. (l27b gives the correction to friction, 

St] = -A'/(0+) [ tRx=o,t - A'^iO^) [ tRx,tg{x) 



A'/(0+)/i+Ag(0+)/2 



(28) 



where we have introduced Ai{u) = Ai(u)/c^. The one-loop 
integrals Ii and I2 diverge logarithmically and for e,5 ^ 






FIG. 3: One-loop diagrams correcting disorder. The dotted line cor- 
responds to either SR disorder vertex (dashed line) or to LR disorder 
vertex (wavy line). Diagrams of type a, b, and c contribute to SR 
disorder. Only diagrams of type d correct the LR disorder. 



read 



i—d 



£ 

,--5 



0{l), (29) 



Q TTl 



where we set m = m/Vc and Kd is the area of a d- 
dimensional sphere divided by (27r)'^. To remove the poles 
in the mobility, we introduce the corresponding Z factor 
?7f{ = Z^^[Ai\rj, which to one-loop order is given by 



1 - A'/(o+)/i - a;,'(o+)/2 



(31) 



In the absence of LR correlated disorder, the elasticity re- 
mains uncorrected to all orders due to the STS, while here the 
correction reads 

5c = ^A'^iO) l^^x'RxM^) 

= -T^A^'(O) / g(g)Vg , ^ , 
2d Jq ^cq^ 

= -c^AUofjm-'. (32) 
d 

We have not set the second derivative at 0+ as A2 remains 
analytic as is discussed below. Furthermore, the correction to 
elasticity (|32] | is finite for e, ^ ^ 0, and thus c does not acquire 
an anomalous dimension. However, we expect corrections at 
two-loop order If this is the case, one has to introduce a Z- 
factor that renormalizes elasticity: cn = Z^^[Ai]c with Zc = 
1 + 0(A|). 

In principle, due to the lack of STS, the KPZ term 
X{\/uxt)'^ breaking the symmetry u —u can be generated 
in the equation of motion (fTSl l at the depinning transition. In- 
deed, diagram e in Fig. |2] when expanding A(u) to second 
order in u, using Eq. ( |27] i, gives 



2d ^ 



(0+) 



x^R 



■x,t. 



,{x). 



(33) 
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Moreover, the term with cubic symmetry (M = 2) and terms 
with higher-order symmetries (M > 2) XM^iidiUxt)'^'^^ 
can be generated by diagram e, 



5\m = 



1 



d(2Af)!^2 



(2A/+1) 



(0+) 



i?..tg(a:)^a;,2*'. (34) 



However, as we will show later, if we start from bare analytic 
disorder distribution, the LR disorder remains analytic along 
the FRG flow and the corresponding FP value [u) is also an 
analytic function. Thus terms ( |33] l and ( |34] l are zero, provided 
that they are absent in the beginning. Moreover, the terms (|34] | 
are UTelevant in the RG sense for M > 2 [but not the KPZ 
term (|33Ti; see Ref. ll46ll 1. This proves that our bare model (ISTTi 
is a minimal model for the description of elastic manifolds in 
a random media with LR correlated disorder. 

The corrections to disorder are given by the diagrams 
shown in Fig. [3] The corresponding expressions read 

5^k^{u) = -{k[{uf + [Ai(u) - Ai(0)]A'/(m)}/i 
-{2k'^{u)K'^{uf + [A2(«) - A2(0)]A'/H 
+Ai(u)A'2'(u)}/2 - [K'^{uf + A2(w)A^'H]/3, (35) 



5^K2{u) = -AiiO)A'^{u)h - A2(0)A'2'(u)J2. 



(36) 



The one-loop integrals Ii and I2 have been defined in 
Eqs. (|29] l and (l30l l. whereas /a is given by 



2(a-d) 



0{l). 



(37) 



,^ (g2 + m2)2 26 -e 
Let us define the renormalized dimensionless disorder A? as 



m-^Af = Ai(w) + <5iAi(w), 
m-'A^^ A2(w) + (5iA2(u). 



(38) 
(39) 



The Pi functions are defined as the derivative of A^(u) with 
respect to the mass m at fixed bare disorder Ai(u). In or- 
der to attain a fixed point, it is necessary to rescale the field 
u by and write the (3 functions for the functions A^ 



2C)Ai(u) + CuA;(w) 



1 d- 



-^[A,{u) + A2{u)f + AA';{u), (40) 
(S - 2C)A2(u) + (uA'^iu) + AA'^iu), (41) 



where A = [Ai(0) + A2(0)] and de ;= -m-^. 

The scaling behavior of the system is controlled by a sta- 
ble fixed point [A^(ii), A2(ii)] of flow equations (ffOl i and 
PTt . To determine the critical exponents, let us start from 
power counting following Ref. |i46il . The quadratic part of ac- 
tion (I2TJ is invariant under x ^ xb, t ^ t¥-, u — > uh'^, 
u ub'^"^-^^^''-^^ . Under this transformation, the mobility, 
elasticity, and disorder scale at the Gaussian FP as c ~ b^ , 

_ b2-z+i,^ _ ^4-d-2C+2V^ and A2 - fe4-a-2C+2V., 

Thus SR disorder becomes relevant for C, — ijj < (4 — rf)/2 
while LR disorder is naively relevant for C, — ifj < {4 — a)/2. 



Note that in the presence of STS, = 0, and we recover the 
conditions obtained at the end of Sec. |III] The actual value 
of C, will be fixed by the disorder correlators at the FP. The 
elasticity exponent ijj and the dynamic exponent z read 



-m— — In Zc{Ai 
dm 



z = 2 + Vj + to— lnZ„(A,) 
am 



(42) 
(43) 



where subscript "0" means a derivative at constant bare pa- 
rameters. To one-loop order this yields 



z = 2- AJ"(0)- A;"(0). 
The scaling relations then read ll46h 



1 



2-C + ^ 



/3 = ^.(z-C) = 



2-C + V' 



(44) 
(45) 



(46) 
(47) 



At zero velocity, the above calculation can be consid- 
ered as a dynamical formulation of the equilibrium problem. 
However, one has to be careful with mapping the dynamic 
FRG equations to the static equations, because as shown in 
Ref. 1 19] the bare relation A^ — —R'l{u) breaks down for 
the SR case at two-loop order. The standard derivation of the 
FRG equations in the statics is based on the renormaUzation 
of the replicated Hamiltonian. We have checked that similar 
to other systems with only SR disorder, the static FRG equa- 
tions for systems with LR disorder can be obtained from the 
dynamic flow equations to one-loop order using the identifi- 
cation Ai — —R'l{u). They read 

diRiiu) = {e~4:C)Ri{u) + CuR[{u) 

+ ^[R'{{u) + R'^iu)]^ + AR'liu), (48) 
diR2{u) = {6~iC)R2{u) + CuR2{u) + AR'^{u), (49) 



R'm- 



where A = -[i?i(0) 

In the case of the model with correlator given by Eq. ( fTsT l. 
one has to distinguish between the transverse and parallel di- 
rections, and therefore introduce corresponding elastic coeffi- 
cients c± and C||. In the transverse direction, disorder is only 
(5-correlated and as a result the transverse elasticity is not cor- 
rected and can be set to 1 . The power counting shows that the 
LR disorder is naively relevant for 61 = 4 — d± — a < 0. 
The one-loop integrals are logarithmically divergent and for 
e, (5i -> are given by Eqs. (|29j, and ^ with 6 ^i. 
Thus the above renormalization can be generalized to model 
( fTSl l if one formally replaces S ^ 61. 

Let us show how a nonanalyticity of the disorder appears in 
the problem. We start from the bare analytic correlators with 
A^'(O) < 0. The flow equation for y -A'/(0) - A^'(O) = 
R'{"{0) + R'^"{0) > Oreads 



diV = £y + iy^ + l{m), 



(50) 
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where 7(m) = (e — 6)^2(0). As we show below, the func- 
tion A2(u) remains analytic along the whole FRG flow and at 
the fixed point (FP). The solution of Eq. (l50l l for any function 
7(771) bounded from below blows up at some finite scale m*, 
which can be associated with the inverse Larkin length. This 
blowup of y corresponds to the generation of a cusp singular- 
ity: Ai(u) becomes nonanalytic at the origin and acquires for 
m < m* a nonzero A[{0^). The precise estimation of the 
Larkin scale requires the solution of the pair of flow equations 
for both Ai{u). 

Before studying different FPs, let us note an important 
property that is vahd under all conditions: if Ai{u) (i = 1, 2) 
is a solution of Eqs. ( |40l l and dTTT i. then K?Ai{u/n) is also 
a solution. Analogously, if Ri{u) is a solution of Eqs. (l48l l 
and ( |49] l, then K^Ri{u/ k) is also a solution. We can use this 
property to fix the amplitude of the function in the nonperi- 
odic case, while for the periodic case the solution is unique as 
the period is fixed. 



VI. NONPERIODIC SYSTEMS: RANDOM BOND 
DISORDER 

In this section, we study the scaling behavior of an elastic 
interface in a disordered environment with LR correlated RB 
disorder. To this aim, we have to find a stationary solution 
(FP) of Eqs. (|48] | and ( |49] l that decays exponentially fast at 
infinity as expected for RB disorder The SR RB FP with 
R2 {u) = 0, which describes systems with only SR correlated 
disorder, was computed numerically in Refs. fl2l[T7llT8ll . The 
corresponding roughness exponent to one-loop order is given 
by CsR = 0.208 298£ + 0(£2). We now look for a LR RB FP 
with R2 (u) 0. Integrating Eq. ( |49] |. we obtain 



di / R2{u) = ((5 -50 / R2{u) 



(51) 



Therefore, the new LR RB FP, if it exists, has 



Clrrb = 7 + 0(e^<5^£<5). 
5 



(52) 



The direct inspection of diagrams contributing to the FRG 
equation for R2 shows that the higher orders can only be lin- 
ear in even derivatives of i?2(w). The only term that is linear 
in R2{u) comes from the renormalization of the elasticity and 
can be rewritten as 2%l)R2 {u) to all orders. Therefore, in higher 
orders we have 



dt / R2{u) = ((5-5C + 2V') / i?2H, (53) 



and as a consequence, duR2 {u) is exactly preserved along 
the FRG resulting in the exact identity 



Clrrb 



5 + 2il} 



Using our freedom to rescale Ri{u), we introduce S := 5/e, 
Ri{u) =: er,{u) and fix r'{{0) =: -x and r^'(O) = -1, 
where x is the parameter to be determined. The stationarity 
condition of Eqs. (l48T l and ( |49] l reads 



4 -\ 5 I 

1 " -^5jri{u) + -ur[{u) 

+ \[r'i{u) + r'i{u)Y + (1 + xYiiu) = 0, (55) 
3 S 

(u) + -ur^(ii) + (1 + x)r'^{u) = 0. (56) 

5 5 



Since Eq. ( |56] | is linear in r2, it can be solved for fixed x by 



r2iu) 



5(1 + a;) 



■ exp 



10(1 



(57) 



From the Taylor expansion of Eq. ( fSSl l around u = 0, we find 

5(1 



ri{0) = 
r[{0) = 0, 
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(58) 



where the second condition excludes the possibility of a super- 
cusp (the first line does not diverge since a; = 1 for S ~ 5/4). 
Thus for fixed 6 the simultaneous equations ( fSST l and ( |56] ) have 
a unique solution for any x, but only for a specific x does the 
solution ri (u) decay exponentially fast to for large u. To de- 
termine this value, we employ the shooting method choosing 
X as our shooting parameter For fixed x, we integrate numer- 
ically Eq. (|55] l with r2{u) given by Eq. dSTT i from to some 
large Umax with initial conditions (|58] l. Then the shooting pa- 
rameter X can be found by solving numerically the algebraic 
equation ri{u^ax]x) — 0. Increasing Umax, we acquire the 
desired accuracy for x and ri (x). We were able to find the nu- 
merical solution with reasonable accuracy only for S > 1.1. 
The typical FP functions r J (u) and (u) are shown in Fig. |4] 
The actual values of x obtained by shooting for different S are 
summarized in Table H] 

Let us now check the stability of SR and LR FPs. To 
that end we linearize the FRG equations about each FP. In 
the vicinity of a FP, the linearized flow equations have solu- 
tions that are pure power laws in to, i.e., scale as m^^ with a 
discrete spectrum of eigenvalues A. A stable fixed point has 
all eigenvalues A < 0. Substituting r.^ r*{u) + Zi{u) 
into the flow equations and keeping only terms that are lin- 
ear in Zi{u), we derive the linearized flow equations at the FP 
{rliu),r*2{u)}, 

(1 - AC,)z,iu) + Ciuz[{u) + [rr'iu)+rnu)] 
X [zUu) + z!^{u)] + {l + x)z'{{u) 
+ Aorl"{u)^Xzi{u), (59) 

[S - 4Ci)z2(u) + Ciu4(w) 

+ (1 + x)z'^iu) + Aor*2"{u) = Xz2{u), (60) 



(54) 



where we have introduced ( — eCi, Aq 
and A is also measured in units of e 



-K'(o) + 4'(o)], 

Note that because 
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TABLE I: Long-range correlated random bond fixed point. Shooting 
parameter x = — r"(0), the maximal eigenvalue and the universal 
amplitude for different values of S. 















1.1 


1.931986 




33.89 


1.2 


1.121722 


-0.160 


31.37 


1.3 


0.922046 


-0.262 


31.64 


1.4 


0.825747 


-0.365 


32.41 


1.5" 


0.766976 


-0.469 


33.34 


2.0* 


0.639151 


-1 


38.44 


3.0"^ 


0.562357 


-2.120 


48.10 


oo 


0.463619 




oo 



"Random lines in a planar interface (d = 2, a = 1). 
'Random lines in a 3d manifold (d = 3, a = 2). 
■^Random planes in a 3d manifold {d = 3, a = 1). 



of the freedom to rescale ri{u), we always have the eigen- 

mode zp"* with marginal eigenvalue Aq = 0. As shown in 
Ref. I.45l1 for SR RB FP the corresponding eigenfunction is 



given by zf^ = ur*'(u) — Ar'l{u), z'^' — 0, while the next 

0} is the SR 



(0) 



(1) _ 



eigenvalue Ai = — 1 corresponds to z\ 

(1 - 4Cp)rt(u),z^°^ = 0. Here {r^r^ 
RB FP and the Taylor expansion of the function rj' can be 
found in Ref. |45|]. Thus the SR RB FP is stable in the 
SR disorder subspace (r2 = Z2 — 0). Let us check its 
stability with respect to introduction of LR correlated disor- 
der. From Eq. ( l60l l it follows that the maximal eigenvalue 
Amax = 5 — 5^1 ^ corresponds to the exponential eigenfunc- 
tion Z2{u) = exp(-CpwV|2rt"(0)|) withrt"(0) = -0.577 
for SR RB FP. As a consequence, the LR correlated disorder 
destabilizes the SR RB FP if (5 > 5Cp w 1.041, or equiva- 
lently, using Eq. ( |52] |. if C,^^ < . This criterion was of 
course expected. 

We now check the stability of the LR RB FP 
{r\(u),r2{u) ^ 0}. It also has a marginal eigenvalue Aq = 



with eigenfunctions given by z. 



(0) _ 



ur*'{u) - ^r*{u) that 



can be checked by direct substitution into Eqs. ( |59] l and ( |60] |. 
Equation ( l60l l allows for an analytical solution that reads 



Z2{u) 



5Ao 
2<5 + 5A 



5(1 



exp 



10(1 + a;) 



(61) 



We are free to fix the length of the eigenvectors, for instance, 
by the condition 2:2(0) = 1, which gives 



Ao = ^{l + x)[25 + b\]. 
3(5 



(62) 



Thus to find the eigenvalue A and the eigenfunction zi, we 
have to solve Eq. ( |59] ) with condition z"(0) = —1 — ^0 and 
require an exponentially fast decay of the solution at large u. 
The only case for which we succeeded to construct the solu- 
tion analytically \s 5 — 2, which is depicted in Fig.|5] It has 




6- - ■ 
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FIG. 4: Fixed point describing the interface in a medium with long- 
range correlated random bond disorder (LR RB FP) for 5 = 2. The 
SR part rj" (u) is a nonanalytic function with rl"{Q^) 7^ 0. The LR 
part rj (u) is an analytic function. Here we report minus their second 
derivative. 



eigenvalue A = — 1 and reads 

Zl{u) : 
Z2{u) 



^_url'{u) + Kl{u) + Kl{u), (63) 



1 



ur*^ {u) + r*2{u)]. 



(64) 



For other values of 5 we solve Eq. ( |59l ) numerically using A as 
a shooting parameter and require an exponentially fast decay 
of zi {u) for large u. To compute the numerical solution, we 
need the initial conditions. Expanding Eq. ( |59l ) in a Taylor 
series, we obtain 



zi(0) 

z'M 



5[a;2(2(5 + 5A) + bx{5 + A) + 3(5] 



3(5(5 - 4(5 - 5A) 



0. 



(65) 
(66) 



Apart from the marginal eigenvalue Aq = 0, the largest eigen- 
value is Ai. It is shown for different 5 > 1.1 in Table U The 
negative sign of Ai reflects the stability of the LR RB FP. For 
S < 1.1 we failed to compute the numerical solution with rea- 
sonable accuracy. However, the largest eigenvalue computed 
at LR RB FP Ai tends to for ,5 ^ 1.1 and the SR RB FP be- 
comes unstable for S > 1.041 with respect to LR-correlated 
disorder Thus we expect that the LR RB FP is stable for 
S > 1.041. Moreover, the largest eigenvalue within accessible 
accuracy is well approximated by Ai = 0.1917(Ci^ — Ci"^ )' 
which gives Ai = —0.06 for (5 — 1.1. 

Besides the roughness exponent, there is another universal 
quantity that is of interest. This is the displacement correlation 
function, which behaves like 



-{d+2C) 



(67) 



Let us show that in contrast to systems with only SR- 
correlated disorder, this system, whose behavior is controlled 
by the LR RB FP, has a universal amplitude Ad- Indeed, 
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-0.4- 



-0.8 



-1 7- 

1 ' 



^7 



FIG. 5: Second derivative of eigenfunctions 21 (it) and 22 (f) com- 
puted at the LR RB FP for 5 = 2. 



according to Eq. ( fSTT ). the integral J duR2{u) is preserved 
along the flow and is fixed to its bare value Q, where we 
have introduced the actual renormalized correlator i?2 which 
is connected to R2 given by Eq. ( |57l ) by the relation i?f = 
K'^Ri{u/ k) with k given by 



)l/5 



3/10 



(27r)i/io U(l + x) 



(68) 



where we used / duB^{u) = Q. Then the amplitude can be 
written to one-loop order as follows lfl9ll : 

A4 



A4 



where we have introduced the universal function 



3/5 



(69) 



(70) 



Values for x(S) and for B{S) for different 5 are shown in Ta- 
ble |D 



VII. NONPERIODIC SYSTEMS: RANDOM FIELD 
DISORDER 

We now address the problem of an elastic interface in a 
medium with LR-correlated RF disorder. We expect that sim- 
ilar to systems with uncorrected disorder, this universality 
class also describes the depinning transition. To see that sys- 
tems with RB disorder flow in the dynamics to the RF FP, one 
has to include either effects of a finite velocity or consider 
two-loop contributions, which go beyond of the scope of the 
present work; but we expect the mechanism to be the same as 
inRef. 1,18.1 . 



Let us look for a solution of Eqs. ( |40] | and ( l4TT l. which de- 
cays exponentially fast at infinity as expected for RF disorder. 
From Eq. (|4TI ) it follows that (hereafter we drop the tilde on 

A,) 



di J A2iu) - [6-30 J A2H. (71) 


Therefore, A2(u) is preserved along the FRG flow fixing 
the roughness exponent to 



Clrrf- 3+0(£',^',£'5), 



(72) 



which coincides with the Flory estimate. Introducing 
Ai{u) = ey.,{u), C = eCi and fixing yi(0) = x, ^2(0) = 1, 
we can rewrite the stationary form of Eqs. ( |40] i and i4T[ as 
follows (Ci = S/3): 

1 (f 

(l-2Ci)yiH + Ciuy[{u) - --^[yi{u) + y2{u)]^ 

+ [l + x]y';{u)=0, (73) 
(5-2Ci)y2H + CiU2/2(") + [l + ^]2/2(«) = 0. (74) 



Equation ( |74] | can be solved analytically giving 

2/2 (u) = exp 



6(1 + x) 



(75) 



Substituting the FP function (iTSl l in Eq. ( f73] l. we obtain a 
closed differential equation for the function yi{u). Expand- 
ing around w = 0, we find 



yi(o) 
yi(o) 

^2(0) 
^2(0) 



1 



9x + 3S — 6xS, 
2 



3 

1 Sx 
3 ~ 9x + l' 
0, 



3{x+l)' 



(76) 

(77) 
(78) 

(79) 



Thus we can compute numerically the solution yi (u) for any 
fixed 6 and j/i (0) = x, however only for special values of x 
does this solution decay exponentially at infinity. The corre- 
sponding solution can be computed using the shooting method 
as described above, using a; as a shooting parameter (see Ta- 
bleHB. 

A pair of typical FP functions is shown in Fig. |6] Sur- 
prisingly, the function yi{u) obtained by shooting satisfies 
du yi{u) — 0, characteristic for RB-type correlations 
along the u direction. In other words, the LR RF FP is in 
fact of mixed type: RB for the SR part and RF for the LR part 
of the disorder correlator. This can be understood as follows: 
Consider the flow of du yi {u). It is obtained by integrat- 
ing the l.h.s. of Eq. (7% from to infinity, and by inserting 
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TABLE n: Long-range correlated random field disorder. Shooting 
parameter x = t/i(0), the maximal eigenvalue and the universal am- 
plitude for different values of 5. 



S x(5) Ai B{S) 

1.1 0.562872 -0.1 140.43 

1.2 0.525082 -0.2 142.23 

1.3 0.496948 -0.3 144.27 

1.4 0.475110 -0.4 146.44 
1.5° 0.457638 -0.5 148.66 
2.0*" 0.404989 -1.0 159.66 
3.0' 0.362329 -2.0 179.04 



"Random lines in planar interface (d = 2, a = 1). 
''Random lines in a 3d manifold (d = 3, a = 2). 
■^Random planes in a 3d manifold (d = 3, a = 1). 



Ci = S/3: 




FIG. 6: Fixed point describing the interface in a medium with long- 
range correlated random field disorder (LR RF FP) for S — 2. The 
SR correlator yl (u) has a cusp at origin and formally corresponds to 
RB type of correlation in direction u. The LR correlator 1/2 (u) is an 
analytic function of RF type. 



de J duyi{u) = 1^1 - Sj J duyi{u) 

^ [yi{u) + y2{u)f\„-(l + x)y[{Q), (80) 



2du 



where we have used the fact that most terms in the FRG- 
equation (|73] | are total derivatives. Finally, we remark that the 
second line of Eq. dSOl l cancels exactly provided that j/2(u) is 
an analytical function, leaving us with 



d,j^ duyi{u) = [l-6j duyi{u). (81) 



This means that for LR-correlated disorder, i.e., 5 > 1, the 
integral of yi indeed scales to 0. A nontrivial fixed point 
is possible at two-loop order for depinning. We remind the 
reader that in ifTsIl it was shown that at two-loop order and for 
SR-correlated disorder, new terms arise in the ERG equation, 
which do not integrate to 0. Indeed, this is the mechanism 
that leads to a breakdown of the result Csrrf — e/3 at de- 
pinning. The same terms will appear here. We expect that the 
additional diagrams due to LR coiTelations do not exactly can- 
cel these terms, especially since these terms are proportional 
to the derivative at the cusp, and LR disorder will probably 
remain analytic, thus not contribute to the anomalous terms. 
These considerations let us expect that at two-loop order the 
integral of yi (u) will be small, but nonzero. 

Let us finally check the stability of the SR RE EP and 
new LR RE ER At the SR RE EP, the roughness is given by 
Csrrf = e/3, and thus we expect the crossover from the SR 
universality to LR at (5 > e, which follows from the condition 
Csrrf = Clrrf- To check the stability of the EPs, we fol- 
low the strategy used for the RB case and linearize the flow 



equations about the RE EPs. We obtain 

d^ r 

(1 - 2Ci)zi{u) + Ciuz[{u) ~ -^^[yl{u) +y2iu)] 

X [zi{u) + Z2{u)]j + {l + x)z'{{u) 

+ Aoyl"{u)^ \zi{u), (82) 

{S - 2Ci)z2{u)+ClUz'2{u) + {l+x)z2{u) 

+ A„y*2"{u) = Xz2{u), (83) 

where we have introduced Aq = 21 (0) + 2:2(0). Eirst, we 
prove our conclusion on the stability of SR RE EP with respect 
to LR-correlated disorder. To that end, we solve Eq. ( l60l l as- 
suming that Ci = Cf^ = 1/3 and a; = yf^{0) = 2/9. We 
obtain that 22 — exp[— Ci^w^/(2j/J^^)] and the correspond- 
ing eigenvalue Amax — 6—3Cf^. Therefore, indeed the SRRE 
EP becomes unstable with respect to LR disorder for 6 > e. 

We now focus on the stability of the LR RE EP. Analy- 
sis of the linearized flow equations ( |82] | and ( |83] | shows that 
there is at least one eigenvector zf ^ — uA*{u) — 2 A* with 
marginal eigenvalue Ao = 0, which corresponds to the free- 
dom of rescaling. Eor arbitrary A, Eq. (|83] | can be solved ana- 
lytically. 



_ Aq [fa^-3(l + a;)] [ 5u^ 
^'^"^ " 3(2^ + 3A)(l + x)2 '^''P l~6(l + a;) 



(84) 



We are free to fix the length of the eigenvectors, for instance 
by the condition 22(0) = 1, which gives 



Ao = -l{25 + 3X){l + x). 
6 



(85) 



Thus to find the eigenvalue A and the eigenfunction zi, we 
have to solve Eq. ( l82b with condition 21 (0) = ^0 — 1 and 
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require an exponentially fast decay for large u. We need the 
initial conditions that can be found by expanding Eq. ( [82] l in a 
Taylor series. 



zi(0) 



-{26 + 3X){l + x), 



(86) 



z[{0) = -|[9A2 + A(12(5-9)](a; + l) + 52(42; + 7) 

- 5{6x + 9) I /26\/9x-6xd + 3d. (87) 

j = 2 is the only case in which we succeeded to find a com- 
pletely analytical solution (see Fig. |7]i. It has Ai = — 1 and 
reads 



uyr{n)~lyKu)-ly*2{u), (88) 



zi{u) 

Z2{u) = uy2{u) +y2{u). 



(89) 



For other values of 6, Eq. (|89i l remains correct, while to obtain 
zi{u), we solve Eq. ( |82] | numerically, using A as a shooting 
parameter. As can be seen from Table [III the largest eigen- 
value satisfies Xi{S) = 1-6 = 3(Ci^ - Ci"^)- This result can 
be obtained analytically as follows. Integrating Eq. ( |82] | from 
to oo, we get the condition 



f 

Jo 



du zi{u) I 1 — (5 — A 



0, 



(90) 



proving that as long as J zi{u) ^ 0, one has A = 1 — 6. There- 
fore, the LR RF FP is stable for 5 > e. Inserting this value of 
A back into Eq. ( |82] |. we obtain after some simplifications 



Y{u) := y^{u)-y^{0) + y2{u)-y2{Q), 
W{u) := Z2{u)[yi{u) + y2{u)]-yi{0)[zi{Q) + Z2{ 



(91) 



The first equation can be integrated with the result 

W{u)] . 



6 ^ r / \ / \ 

-uzi[u) = — [Y[u)zi[u) 



This is equivalent to 



Zl 



z[Y = -W'. 



(92) 



(93) 



The homogenous equation reads 



[\^z,Y]'=^-^. 
L 1 J 3^ 



Its solution is 



zi{u) 



C 

Yin) 



exp 



ds — 
3 Jo Y{s) 



(94) 



(95) 




FIG. 7: Eigenfunctions zi (u) and Z2 (u) computed at the LR RF FP 
for (5 = 2. 



with some constant C. A solution to the inhomogeneous equa- 
tion is obtained by replacing C by C(w), and inserting the 
latter into Eq. ( |93T l. This yields 



CM 



dtW\t)exp 



ds 



Y{s) 



(96) 



Putting together everything, we obtain 
1 



zi{u) = - 



Yin) Jo 



dtW'{t)exp 



- id 

3 1 'y{s) 



(97) 

Let us now compute the universal amplitude defined by 
Eq. (|67] i. According to Eq. ( ItTI i. the integral J du A2{u) is 
preserved along the ERG flow and can be fixed by its bare 
value Q. The relation between the actual renormalized dis- 
order and the rescaled disorder A2{u) given by Eq. jTSl l 
reads A'"^ = /S.i{u/ n). We obtain 



Q 



1/3 



26 



1/6 



ei/3 l37r(l + a:) 



(98) 



where we have fixed J du/S^2{'^) = Q- Then the amplitude 
can be written as 

1 



A = — [Ar(0) + A^'-(0)] 



en 



+ =£1/3Q2/3^(,5), 



(99) 



with the universal functions given by 



B{5) = Sn\l + xf/^ ( — 



1/3 



(100) 



Values of B{5) for different 6 are shown in Table HIl 

Depinning. We are now in the position to study the depin- 
ning transition, which we expect is controlled by the LR RF 
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TABLE ni: Periodic systems with LR correlated disorder. The shoot- 
ing parameter: yi (0) and two first eigenvalues for different e. 



yi(o) 



Ai 





1/3 

1/2 

2/3° 

0.8 

0.9 



0.00971 
0.01089 
0.01183 
0.01348 
0.01645 
0.02346 



0.0 

0.333 

0.500 

0.667 

0.800 

0.900 



-4.089 
-0.500 
-0.280 
-0.125 
-0.013 



"Corresponds to d 
the plane of a CDW. 



2, a = 1, i.e., line defects (e.g., dislocations) along 



FP. The dynamic critical exponent z defined by Eq. ( |45] | is 
given to one-loop order by 



(101) 



where we have used Eqs. ( fTTl ) and (|79] |. which give y"(0) -I- 
2/2 (0) = 1/3 — (5/9. Other exponents can be computed using 
scaling relations ( l47T i and (l46T l, for example 



l3 = l-^- + ^+0{e\S\eS). 



(102) 



It is remarkable that for S > 3e, the exponent /3 is larger 
than 1, and z is larger than 2. This seems to imply some dif- 
ferent physics - yet to be understood - in the avalanche pro- 
cess, which makes the motion slower near depinning than in 
the SR case. The analyticity of A2 seems to suggest some 
smoother motion at large scale, while short-scale motion re- 
mains jerky and avalanche-like. Finally, note that at the SR RF 
FP, zsR = 2 — 2e/9, /3sr = 1 — e/9, and thus, the exponents 
are continuous functions of e and 6. 



VIII. PERIODIC SYSTEMS 

We now study periodic systems with disorder correlator 
given by Eq. (|4|i, which we can refer to as an XY model 
with LR-correlated defects. The results for CDWs with LR- 
coiTelated disorder defined by correlator ( fT5] l can then be ob- 
tained by substituting 6 ^ Si — 4 — d± — a. It is suffi- 
cient to consider the system with the period fixed to 1 , since 
other systems can be related to the latter using the freedom 
to rescale. As a consequence, the roughness exponent for 
periodic systems is C = 0. At variance with interfaces, we 
introduce reduced parameters according to Ai{u) — Syi{u), 
A — yi{0) + 2/2(0), and e — e/6. Then the fixed-point equa- 
tions can be written as follows: 



2du 



2 [yi{u) + + Ay'{{u) = 0, (103) 



eyi{u) 

y2{u) + Ay'i{u)^Q. (104) 
Equation ( |104| i can be solved analytically. Its solution is 

2/2 = y2(0)cos(2^w), ^ = 1/(2^)2. (105) 



0.020 



0.015 



0.010 



0.005 



0.000 



-0.005 



-0.010 




FIG. 8: (Color online) Fixed point of a periodic system with LR 
correlated disorder. The SR disorder correlator yi (u) computed for 
different values of i. 



Equation ( I103l l can be solved analytically for e — 0, 

yi = yi(0)-Hy2(0)[l-cos(2^u)] 

--^V22/2(0) [l-cos(27ru)]. (106) 

ZTT 

The coefficients yi{Q) are determined by potentiality of the 
Ai, i.e., from the conditions 

1 1 

du yi{u) ^ J du y2{u) = 0, (107) 



and the identity yi(0) + ^2(0) = l/(27r)2. They read 

2/1(0) = l/(27r)2 - 1/64, (108) 
2/2(0) = 1/64. (109) 

For e > 0, Eq. ( I103l l can be written in the following form: 

eyi{u) - \^ Uy,{u) + y2{u)f ~ ^1 = 0, (110) 



where 2/2(u) is given by Eq. ( 1105b with 2/2(0) — l/(27r)^ — 
2/1(0). Expanding Eq. dl 10b in a Taylor series about u = 0, 
we find that 2/i(0) = -^Xl^lnf - 2/i(0)(l - e). Thus for 
any fixed < e < 1 and 2/1 (0) we have only one solution 
yi{u), but only for a specific 2/1 (0) this solution fulfills the 
condition 2/1 (1) = yi(0). To find this value, we employ the 
shooting method using 2/1 (0) as a shooting parameter. The 
values of 2/1 (0) computed for different e are summarized in 
Table |lll] The corresponding eigenfunctions 2/1 (u) are de- 
picted in Fig. [8] 

While the roughness exponent is zero, the system forms a 
Bragg glass phase with a slow growth of the displacements 
according to 



{u^ - Uo)2 = Adl-a \x\, 



(111) 
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where Ad is a universal amplitude, which to one-loop order is 
given by 



A. 



(LR) 



2IQ 



[Alio) 



(112) 



where we have restored the factor of I/K4 previously ab- 
sorbed in Ai . The SR periodic FP is characterized by -A^;^^"* = 
e/18 + ©(e^). It is interesting to compare Eq. ( II 12b with 
the prediction of the Gaussian variational approximation for 

£/(27r^). We expect the 



the SR disorder case A 



(SR) 
d,GVA 



crossover between LR and SR FPs at A 
LR disorder to be relevant for 



(LR) ^ ^(SR) 



I.e., 



e 9 

T < ^ w 0.912. 



(113) 



We now check the stability of the SR and LR periodic FPs. 
The flow equations linearized about the FP read 



r 

sziiu) - -^^[yl{u) +y^{u)][zi{u) + Z2( 

+ Az'{{u)+Aoyl"{u) = Xzi{u), (114) 
Z2{u) + Az'^iu) + Aoyfiu) ^ Xz2{u), (115) 

where we have introduced Aq = zi(0) + 22(0). Let us recall 
that the SR FP is unstable with respect to nonpotential pertur- 
bations even in the subspace of SR disorder. Indeed, the SR 
periodic FP, 



The corresponding eigenvalue Asr = 1 — £7r^/9 confirms our 
estimation for the stability of the SR periodic FP ( II 13l l. For 
the LR periodic FP, we still have Eq. (11 17l l with 



An = - 



X 



l-47r2yi(0)' 



ziiO) = Ao 



1. 



(118) 



Equation ( II 141 ) has a periodic solution only for a discrete set 
of eigenvalues (the first two are shown in Table lllll i. It 
follows from the table that Ai = e > 0. In analogy with the 
SR periodic FP, the LR periodic FP is unstable with respect to 
a nonpotential perturbation corresponding to Ai. The latter is 
obtained by integrating Eq. ( II 151 ) over one period. 



(e — A) / du zi{u) — 0. 



(119) 



As long as the integral does not vanish, this gives the reported 
eigenvalue Ai — e. Indeed, as can be seen from Fig. |9] we 
have duz[^'' {u) ^ and duz[^^^ (u) ^ for n > 2. 

Depinning. We now focus on the depinning transition of 
the periodic system with LR-correlated disorder. At the LR 
periodic FP, we have 



yi(0) = l + --47rV(0), 



(120) 
(121) 



and thus (0) + 1/2 (0) = The dynamic critical exponent 
z defined by Eq. ( 1451 ) reads to one-loop order 



At(u) 



1 

--u[\-u) 



A*(u) = 0, 



(116) 



has in the SR subspace the positive eigenvalue Ai = 1, corre- 
sponding to the nonpotential eigenfunction zi = 1. All other 
eigenfunctions are potential, i.e., fulfill condition (11071 ). and 
have negative eigenvalues If we add LR-correlated dis- 
order, the solution of Eq. dl 151 ) yields 



(117) 



Z2{u) = COS 2lTU. 



,LR 



2-- + 0{e\6\e5) 



Therefore, for periodic systems z 



LR 



(122) 
to one-loop order 



IX. FULLY ISOTROPIC EXTENDED DEFECTS 

In this section, we briefly examine the effect of a defect 
distribution isotropic in the whole (x, u) space. Consider first 
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0.6 I •■ •■ •■ •■ •- 

-0.6 -0.4 -0.2 0.2 0.4 

P 

FIG. 1 1 : The roughness exponent of the optimal paths on the plane 
with isotropically correlated disorder (data taken from Ref. |48]). 
Solid line — 3/(4 — 2p) > ^sr ~ 2/3 is the roughness expo- 
nent of the elastic string on the plane with fully isotropic long-range 
correlated disorder. Dashed line = (3 + 2p)/5 > C^sr is the 
roughness exponent of the elastic string on the plane with disorder 
correlated only along the string. 

interfaces in random bond type disorder From Eqs. (|5]l and 
(|9]l, one finds 



R{x,u) = VnB{x,u)VMO,0) ~ , ^ ^\ia/2 ' (^23) 

and thus the u and x dependences are coupled in the bare cor- 
relator For the present discussion, we consider N, the number 
of components of u, arbitrary, hence D = N + d. We recall 
that a — D — Ed- The question of to which universality class 
this model belongs is subtle. It turns out that it does not corre- 
spond to LR disorder in internal space, but rather SR disorder 
in internal space and LR disorder in the u direction, hence 
i?2 = but Ri (u) long range in u. To see this, let us consider 
at fixed u the integral / d'^xR{x, u). We can distinguish two 
cases, as follows. 

(i) For a > d, this integral is convergent at large x, 
hence we clearly have SR disorder in the x direction, and 
Ri{u) ^ lul"^^" at large u. This, however, is LR disorder 
in u. This case has been studied using FRG and yields, for 
a < add, N), a roughness exponent given by the Flory value 
C(a, d) = (4 - d) /(4 + a - d). The value add, N) can be 
estimated using the value for the SR disorder roughness ex- 
ponent, by requiring ({a^d, N),dy^ CsKRB{d,N) (small 
deviations can arise as discussed in [7]). 

(ii) For a < d, the situation is more subtle and one may 
be tempted to argue, since / d'^xR{x, u) diverges in the in- 
frared, that disorder LR in x is produced. This is, however, 
not the case, as can be seen on the Fourier transform R{q, P), 
where P is the momentum associated with u, and q with 
X. One has R{q,P) - (g^ + p2ya-d-Ar)/2^ .^j^j^j^ j^^^ ^ 
well defined Hmit R{q = 0, P) = pa~d-N ^ -pj^j^ 
sponds again, as we argue, to a SR correlator in space with 
Ri(u) - i?i(0) - \u\'^-°-. As is often the case, LR models 
require some trivial subtractions. The subtracted correlator 



R{x, u) — R{x, 0) has indeed a convergent integral ~ 
at large u, while subtracting a u-independent piece does not 
change the model. The critical case a = d is described by 
the logarithmic model R{x, u) — R(x, 0) ~ In \u\, which has 
( = (4 - d)/4 in all dimensions [47]. 

To summarize, isotropic distributions of defects isotropic in 
the [x, u) space also yield LR models, but not of the type (|4|i 
studied here. For isotropic line defects, one finds ^ = (4 — 
d)/(3+iV) (i.e., C = 3/4 for a directed polymer in D = 1 + 1, 
( = 3/5 in D = 1 + 2, and for an interface D = 2 + I, 
C = 2/5). Isotropic planar defects yield ( = (4-d)/(2 + iV), 
hence ^ = 2/3 for a. [D = 2 + 1) -dimensional interface. 
This case is illustrated in Fig.[TO] Note that in that case there 
are infinitely many lines of defects inside the interface with 
random directions (the intersections of the planar defects with 
the interface gives lines), but that this does not suffice to create 
power-law correlations in internal space, as can be seen from 
the example in which the planar defects are orthogonal to the 
interface. 

In Ref. ll48ll . the universal properties of the optimal paths on 
the plane with isotropically correlated random potential which 
coiTelation decays as r'^P^^ were studied using numerical sim- 
ulations. This model is believed to belong to the same univer- 
sality class as the one-dimensional {d ~ 1) elastic interface 
in a medium with fully isotropic long-range correlated disor- 
der with a = 1 — 2p. The roughness exponent of the shortest 
paths computed in Ref. [48] for different values of p is shown 
in Fig. [TT] In the optimal path model, the elasticity is gen- 
erated by disorder and the long-range correlated disorder can 
generate the long-range elasticity, which would decrease the 
path roughness and account for the deviation from our pre- 
diction. The effective elasticity is expected to be described 
in terms of the exponent ip, which computation requires the 
two-loop consideration. 

Finally in the periodic case, such as for CDWs, isotropic 
disorder in the full space (xy , x±) again leads to correlations 
( fTSl ). but now the function g{x^^ — a;'|) decays exponentially 
beyond a length scale set by the disorder period (as can be seen 
in Fourier space considering the discrete P modes). Hence the 
problem is described by the standard (SR) random periodic 
class. 



X. CONCLUSION 

We have studied elastic interfaces and periodic systems in a 
medium with LR coiTelated disorder, both in equilibrium and 
at the depinning transition. This type of long-range correla- 
tion exists in the internal space of the manifold, and we have 
discussed how it can be realized in terms of extended defects, 
or anisotropic defects with a broad distribution of lengths. Us- 
ing a dynamic formalism, we derived the FRG flow equations 
for the SR and LR parts of the disorder correlator and found 
three new FPs, which describe three new universality classes. 
All new FPs are characterized by a nonanalytic SR part of the 
disorder correlator and an analytic LR part. We have com- 
puted the corresponding exponents and universal amplitudes 
in a double expansion in e — 4 — d and S — 4 — a. For RB 
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type of disorder, we find that the LR correlation of disorder 
is relevant for S > 1.041e and results in the roughness expo- 
nent C, — (5/5, while for S < 1.041e the scaling behavior is 
conti-olled by the SR RB FP with C = 0.208 298e. We find 
that the presence of RF disorder results in a mixed FP with 
the SR correlator corresponding formally to RB type of disor- 
der and an analytic RF LR correlator. The LR RF FP, which 
is also expected to control the depinning transition, is stable 
for (5 > e giving C = 5/3 and /3 = 1 - e/6 + 6/18. The LR 
correlated periodic FP is stable for e < 0.912(5 and gives a 
slow logarithmic growth of displacements with universal am- 
plitude A]}^ = (5^/(27r^). It is remarkable that this type of 
disorder yields an exponent f3 for the velocity-force character- 
istics that can be larger than unity and a dynamical exponent 
larger than 2. This striking behavior might be relevant for ex- 
periments, and gives a strong motivation for numerical studies 
of the problem, e.g., to understand the nature of motion at the 



depinning transition in these systems. 
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